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A way to explain the puzzling difference between the pion form factor as measured in e^e~ 
annihilations and in r decays is discussed. We show that isospin symmetry breaking, beside 
the already identified effects, produces also a full mixing between the p", lu and cj) mesons 
which generates an isospin component inside the p" meson. This effect, not accounted for in 
current treatments of the problem, seems able to account for the apparent mismatch between 
e^e~ and r data below the S mass. 



1 Introduction 



In order to get a theoretical estimate of the muon anomalous magnetic moment g — 2, one 
needs to estimate precisely the photon vacuum polarization (see Jegerlehner^ for a comprehensive 
review). Its leptonic part can be computed theoretically to a high precision from QED, but the 
dominance of non-pertubative effects in the low energy region prevents to perform likewise 
starting from QCD in order to estimate the hadronic part. This is instead done by means of 
a dispersion integral involving the measured cross section a{e~^e~ — > hadrons); however, the 
integration kernel is such that the low energy region contribution is enhanced by a ~ 
factor. Because of this, the non-pertubative region provides, by far, the largest contribution to 
the hadronic vacuum polarization (VP). Additionally, the annihilation process e+e~ — > tt+tt" 
alone happens to provide more than 60 % of the total hadronic VP. 

As one has to rely on data in order to estimate this contribution, the precision of the 
measured a{e~^e~ tt+tt") is clearly an important issue. For this purpose, several sets of data 
collected in e^e~ annihilations at low energies during a long period of time are available. This 
covers the former data sets collected by the OLYA, CMD and DM1 Collaborations - which are 



gathered in the review by Barkov et al'^ and the data sets more recently collected by the 
CMD2^'^'^ and SND^ Collaborations. Additional data sets taking advantage of the initial state 
radiation mechanism have also been collected by the KLOE, BaBar and Belle Collaborations 
and are expected to become available soon. 

Moreover, high statistics data on the decay u-rn^n^ are also available from the 

ALEPH ^ and CLEO Collaborations. As the pion form factor [i.e. the dipion mass spec- 
trum) in r decays and in e^e~ annilhilations are related by the Conserved Vector Current 
(CVC) assumption, these data are expected to be useful in order to improve the estimate of the 
photon hadronic VP. Indeed, these two kinds of data can only differ by isospin breaking effects 
which are subject to accurate estimates. 

Isospin symmetry breaking effects have been especially studied in order to include the r 
data in the estimation of the photon hadronic VP. This covers non trivial effects specific of the r 
decay like the short range ^ and long range ^'^'^^ isospin breaking factors, but also more standard 
effects easier to account for : mass differences between charged and neutral pions and charged 
and neutral p mesons, the — p^ width difference, or the to and (p contributions to the e~^e~ 
annihilation amplitude (see, for instance, the work by Davier et al. ^^'^^). 

As a preliminary step in the process of including r data in estimating the photon hadronic 
VP, the comparision has been done of the pion form factor as measured in e~^e~ annihilations 
and as derived from r decays while accounting for all known isospin breaking effects appro- 
priately I2,i>i3 'pjjjg comparison, however, clearly exhibits an unexpected s-dependence of the 
difference between the e+e" data and the pion form factor function reconstructed from r data, 
as reported still recently by Davier °. 

This mismatch is an important issue as the photon hadronic VP as reconstructed from 
e+e~ data leads to a theoretical prediction for the muon y — 2 at ~ 3.3 cr from its measured 
value in contrast, the prediction derived from r data fed with all currently identified isospin 
symmetry breaking effects provides an expectation in close agreement " with the g — 2 value 
directly measured at BNL Therefore, the question is whether the {e'^e~ — r) mismatch can 
be explained by physics only connected with isospin symmetry breaking — and then something is 
missing - or if it calls for another kind of physics effect (actually hard to identify) . Responding 
this question by leaving e~^e~ data beyond any doubt may point towards a new physics effect 
exhibited by the muon anomalous magnetic moment. 

2 A missing piece of isospin symmetry breaking ? 

As it is clear that all identified (and listed above) effects produced by isospin symmetry 
breaking should be considered, the question is about a possible missing piece in the isospin 
symmetry breaking procedure (or a piece not appropriately accounted for). 

Actually, a clue towards the solution we propose has been given by Maltman using sum 
rules derived from an OPE input, he concluded that the p part of the e'^e~ form factor data 
was inconsistent with being isospin 1, in contrast with the corresponding information provided 
by r data. This statement implies that either the quality of the available e'^e~ data** can be 
questioned or that the p^ meson is not a (pure) isospin 1 object. 

Up to now, the model amplitudes used to describe the neutral and charged p mesons may 
differ by feeding their propagators with different masses and widths to be fit with data; of 
course, the pion mass difference is also fed in, together with the u (and <p) meson(s) propagators, 
generally Brcit Wigner formulae. However, as the e~^e~ data clearly exhibit ^'^'^'^ the narrow 
(isospin 1 part of the) uj interfering with the broad p^, one may ask oneself about the existence 
of a (broad) isospin part of the p^ meson which might make it differing from its charged 



"See also A. Hoecker, http://moriond.in2p3.fr/QCD/2008/MorQCD08Prog.html, Moriond QCD, March 2008. 
''One may note that these data has now been confirmed by several new and precise measurements. 



partner beyond genuine mass effects. Stated otherwise, the question is wlicthcr mass and width 
differences for the p mesons exhaust isospin symmetry breaking in the pion form factor. 



3 The Pion Form Factor at One Loop 

In order to make our statements exphcit, we have found it appropriate to work in the 
framework of the Hidden Local Symmetry (HLS) model ^^'^^ However, our arguments are 
certainly valid in most other VMD-like models. In the HLS model, the pion form factor for 
both e^e~ annihilation and r decay writes : 

i^.(^) = (1 - |) - {fp = a9fl , W = f) , (1) 

where a is a parameter specific of this model - close to 2 - and g is the universal vector cou- 
pling - close to 5.5 (from QCD sum rules). Dv{s) = s — mj, is the inverse p bare propagator 
(m-p = ag"^/^). While including one loop effects, Dy^s) acquires a pion (and kaon) loop term 
Iip{s) which shifts the p pole off the real s-axis. The transition amplitude from 'y/W to (neu- 
tral/charged) p is also dressed by loop effects; this turns out to perform the change : 

fp^Fp = f,-Uw/^is) (2) 

in the expression'^ for F.„-{s). The 3 loop functions ^p{s) and n^/^(s) just defined fulfill each 
a dispersion relation and their imaginary parts are influenced by SU(3) flavor symmetry 
breaking. Each of these carries a subtraction polynomial, which has been chosen of degree 2 and 
vanishing at the origin. Additionally, it has been possible to relate the subtraction polynomials 
for Ilw{s) and n-y(s). All this lessens significantly the model parameter freedom. 

Breaking isospin symmetry turns out to multiply [^^^^(s)^ by some specific factors ^'^'^'^"'^ not 
discussed here (see and add the co and (p contributions to F^{s). This, however, has been 
shown insufficient in order to restore consistency between e+e~ and r data-*^^'^^'-*^^. 



4 A model for p^, ui, (j) mixing 



As it is clear that the isospin symmetry breaking effects listed above have to be taken into 

account, the question is rather about a missing piece in the scheme outlined in Section 3. While 
working at one loop order, the HLS model provides self-masses already referred to for the p 
meson propagators. However, it also contains the piece : 



lag 



(p? +UJI- V2zv(f>i) K-dK+ + (p? -LOI + V2zv(t>i) K'^dK^ 



(3) 



which - through kaon loops - generates transitions'^ among the so-called ideal (bare) fields pj, 
ujj and (pi with no counter part affecting the p^ field. We have : 



^u,<t>{s) = -9ioKKgct,KK[^±{s) +no(s)] 

^pu){s) = 9pKK9u>kk[^±{s) - no(s)] 
^p4>{s) = -9pKKg4,KK[^±{s) - no(s)] 



(4) 



■^In this paper, we only outline the method and refer the reader to ^® for detailed information. 

'^Il-y{s) refers to the pion form factor in e'^e~ annihilation, which will be denoted F^{s), while Ilw{s) refers 
to the pion form factor in r decay correspondingly denoted F^{s). 

*Thc anomalous and the Yang-Mills pieces of the full HLS Lagrangian contribute also to the mechanism we 
outline here in a quite analogous manner these contributions will be skipped from now on. The constants za 
and zv in Eq. 3 are SU(3) breaking parameters which have to be determined by fit. 



as transition amplitudes between the (ideal) p^, oji and (/>/. n-i-(s) and no(.s) denote, resp. the 
charged and neutral kaon loops amputated of the coupling constants factored out for sake of 
clarity. These loops are defined by dispersion integrals over their imaginary parts and contain 
subtraction polynomials {P±{s) and Pq{s)) real for real s, the invariant mass squared flowing 
through the vector meson lines. These polynomials arc chosen of degree 2 and vanishing at 
s = 0; their coefficients have to be fixed by external conditions. If isospin symmetry is conserved 
one may assume that P±{s) = Po{s) and, then, Hf^^s) and 11^0(5) identically vanish; when 
isospin symmetry is broken this condition is certainly no longer fulfilled. Therefore, the HLS 
model which always predicts ujj — 4>i transitions (as Iiuiii){s) never vanishes identically), predicts 
additionally pj — loj and — (t>i transitions when isospin symmetry is broken. 

Therefore, in the general case of isospin symmetry breaking, there are transitions among 
the ideal vector fields. If one defines the physical vector fields as eigenstates of the vector mass 
matrix, as the amplitudes in Eqs. 4 provide non-vanishing entries in the vector meson squared 
mass matrix, these cannot coincide with their ideal partners. Let us define the vector V and V/ 
as the vectors constructed with (resp.) the (physical) p", uj and 4> fields on the one hand, and 
the (ideal) p^) ^"^^ 'Pi fields the other hand. Then the mass eigenstates of the vector meson 
squared mass matrix and the ideal fields are related by F = R{s) Vj and Vj = R{s) V with : 

f ei Mfi 



R 



'^■wk{s) - £2 (1 - zv)^? + (1 - lJ?)e2 

pei pe2 



1 



\ (1 — zy)m2 + n7r7r(s) — /x2e2 {I — zv)m? + {I — p?')e2 

(5) 

where ei = Xlp^{s) and 62 = ^p(f>{s) are functions of s, real below s ~ 1 GeV^. Indeed, the loop 
imaginary parts start at the corresponding two-kaon thresholds. One neglects terms of second 
order in ei and/or €2- n^7r(s) is the pion loop representing the bulk of the p self-energy and 
= ag^f^ is the unperturbed p meson mass squared. 

Performing the change to physical fields into the HLS Lagrangian generates isospin sym- 
metry violating couplings of the lo and (j) fields to vr+vr^, while leaving the p^ coupling to vr+vr^ 
identical to that of its ideal partner at leading (first) order in the e-j. In contrast, the ^ — p^ 
transition amplitude (named fp in Eq. 1) is modified to^^ fp + ^fp{s) at leading order in the e^, 
while the W — p^ transition is obviously unaffected. We thus get, using obvious notations : 



fp = agfl 

n =agfl 



3 n7r7r(s) -62 3 (1 - zv)nn? + TI^t,{s) - p?'e2 



(6) 



Therefore, because of one-loop effects, isospin symmetry breaking introduces a s-depcndent 
difference between the j — p^ and W — p^ transitions; this is entirely due to the fact that ideal 
neutral vector fields cease to coincide with physical neutral vector fields, when defined as mass 
matrix eigenstates. Loop effects always affect the {to, cf)) sector, but the whole (p°, u, (f)) is 
affected only when, additionally, isospin symmetry is broken. Clearly, this effect has not been 
accounted for in previous analyses of the pion form factor in e^e~ and r data. 



5 How to work out the model ? 

The issue now is whether the (p*^, (/>) mixing we just sketched is able to account nu- 
merically for the long standing mismatch between e^e~ and r data. From the point of view of 
data analysis, the number of parameters (coupling constants, U(3)/SU(3) breaking parameters. 



subtraction parameters from dispersion integrals. . . ) in our HLS based model is too large to 
hope fixing them reasonably well using only the e~^e~ and r data. Fortunately, there is a way 
out. 

It indeed happens that the radiative decays {PV'y and P77), which are accounted for by the 
anomalous sector of the HLS Lagrangian, depend on a large part of the parameters involved in 
our model and can serve to fix them quite reliably, even by fitting them in isolation If one 
adds to this data set the leptonic decay information for the to and (p mesons on the one hand, 
and two— pion decay information of the ^ meson •'^ on the other hand, the minimization program 
becomes numerically well defined. This additional data set will be referred to as "decay data". 

Therefore, the resolution method we propose is to consider the e+e" and r data together 
with the decay data. One should stress that^ the form factor F^(s) is entirely determined, from 
a numerical point of view, by the e+e" and decay data in isolation, since actually all parameters 
it depends on are already involved in the decay widths considered or in -F^(s). Stated otherwise, 
F^{s) can be predicted from our model using only the e+e~ and decay data. We actually consider 
this last property as the main test of validity of our approach. 

Before closing this Section, one should mention that the HLS model, as currently known, 
involves the pseudoscalar meson (P) and only the lowest lying vector meson nonet {V). This 
prevents, for the time being, to include in fit procedures invariant mass region influenced by 
higher mass vector mesons, such as the p(1450) or p(1700). Therefore, the fits to form factors 
will be restricted to the s < 1 GeV^ region, where the corresponding effects should be limited. 

6 Brief analysis of fit results 

Detailed fit information can be found in Ref. where they are lengthily presented and 
discussed. Here, we limit ourselves to the most relevant. One should also mention that data on 
the pion form factor in the close spacelike region ^^'^^ are included in our fits. 

Table 1: Global fit information with special emphasis on pion form factor data subset contributions. 





Full Fit 


Excluding r data 


No — mass shift 


x'/dof 
Probability 


313.83/331 
74.4% 


257.73/274 
75.2% 


321.75/332 
64.7% 






All TimcHkc Data 


187.15/(209) 


170.70/(209) 


192.38/(209) 


T ALEPH 
(x^/points) 


23.86/(33) 


42.27/(33) 


24.28/(33) 


r CLEO 
(xVpoints) 


26.06/(25) 


26.16/(25) 


28.55/(25) 



Table 1 clearly shows that the description of the global data set is quite satisfactory. 
The second data column gives mostly the distance of the model to the r data points left out 
from the fit procedure; this clearly illustrates that F^(s) is indeed numerically derived from the 
HLS model together with data independent of the r form factor. One can even remark that 
CLEO data^ are as well accounted for as when including them in the fit procedure ! Comparing 
the first and tliird data columns, one may also remark that, allowing a different mass for the 

^Thc e+e^ spectrum in the region of the meson is not available as such in the data published by the various 
Collaborations at Novosibirsk. 

^Except for a parameter Sruy which may accout for a (possible) mass difference between p° and p*. 



charged and neutral p mesons provides a marginal improvement. Actually, the corresponding 
mass difference (at the edge of statistical significance) is visible only in ALEPH^ data and needs 
confirmation by forthcoming data sets. 

Therefore, one may conclude that introducing the effects of isospin symmetry breaking (a 
non— zero ei(s), substantially) on vector meson mixing, together with the already reported effects, 
is enough to reconcile the e+e^ and r data. A missing piece in the current isospin symmetry 
breaking procedure is then identified as the effects of the isospin component of the meson 
which has no counter part inside the meson. 



o) Aleph 





V(s) (GeV) V(s) (GeV) 

Figure 1: Distributions of fit residuals for the ALEPH and CLEO t data from the full fit with an allowed mass 
difference between charged and neutral p mesons. The arrows indicate the upper limit of the fit region. The insets 

magnify the p peak invariant mass region. 

Other information is provided by Figs. 1 and 2 which exhibit the fit residuals. One can 

clearly consider them as structureless in the region below 0.9 GeV as no obvious s-dependent 
effect can be observed. One also clearly sees the effects of higher mass vector mesons starting 
as early as around the GeV region. 




V(s) (GeV) 

Figure 2: Residual distribution for all the e^e~ new timelike data over the whole invariant mass interval. The 

inset magnifies the p peak invariant mass region. 



As final conclusion, one may indeed consider that e^e and r data do not exhibit any mis- 



match once all consequences of isospin symmetry breaking arc indeed considered, including the 
isospin component generated inside the meson. Then, it follows from this work that the 
predicted value of the muon anomalous moment derived using e+e^ data is indeed reliable and 
that the actual mismatch is between the prediction of the muon g — 2 and its direct (BNL) 
measurement rather than between e^e" and r data. Therefore, getting an improved mea- 
surement of the muon anomalous magnetic moment becomes a key issue, possibly a window 
on some New Physics. 
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